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Abstract. This paper deals with global convergence to equilibria, and in particular 
Hirsch's generic convergence theorem for strongly monotone systems, for singular 
04 ' perturbations of monotone systems. 
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(T^ ■ 1 Introduction 

^. 

^^ • Monotone systems constitute a rich class of models for which global and 

almost-global convergence properties can be established. They are particu- 
larly useful in biochemical models (see discussion and references in [14, 15]), 
1-^ I and also appear in areas like coordination ([11]) and other problems in control 

C^ • theory ([1]). This paper studies extensions, using geometric singular pertur- 

C I bation theory, of Hirsch's generic convergence theorem for monotone systems 

([4, 5, 6, 13]). Informally stated, Hirsch's result says that almost every bounded 

solution of a strongly monotone system converges to the set of equilibria. 

There is a rich literature regarding the application of this powerful theorem, 

i-j I as well as of other results dealing with everywhere convergence when equilibria 

are unique ([13, 2, 7]), to models of biochemical systems. Unfortunately, many 
models in biology are not monotone. In order to address this drawback (as well 
as to study properties of large systems which are monotone but which are hard 
to analyze in their entirety), a recent line of work introduced an input/output 
approach that is based on the analysis of interconnections of monotone sys- 
tems. For example, the approach allows one to view a non-monotone system as 
a "negative" feedback loop of monotone open-loop systems, thus leading to re- 
sults on global stability and the emergence of oscillations under transmission 
delays, and to the construction of relaxation oscillators by slow adaptation 
rules on feedback gains. See [14, 15] for expositions and many references. The 
present paper is in the same character. 

Our motivation arose from the observation that time-scale separation may 
also lead to monotonicity. This point of view is of special interest in the 
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context of biochemical systems; for example, Michaelis Menten kinetics are 
mathematically justified as singularly perturbed versions of mass action ki- 
netics. A system that is not monotone may become monotone once that fast 
variables are replaced by their steady-state values. A trivial linear example 
that illustrates this point is x — —x—y, ey——y+x, with e>0. This system 
is not monotone with respect to any orthant cone. On the other hand, for 
e <C 1, the fast variable y tracks x, so the slow dynamics is well-approximated 
by i; = ~2x (which is strongly monotone, because every scalar system is). 

We consider systems x = f{x,y), ey = g{x,y) for which the reduced 
system x = f{x, h{x)) is strongly monotone (in fact, a slightly stronger tech- 
nical condition on derivatives is assumed) and the fast system y — g{x,y) 
has a unique globally asymptotically stable steady state y = h{x) for each 
X, and satisfies an input to state stability type of property with respect to x. 
One may expect that the original system inherits global (generic) convergence 
properties, at least for all £>0 small enough, and this is indeed the object of 
our study. This question may be approached in several ways. One may view 
y — h{x) as an input to the slow system, and appeal to the theory of asymp- 
totically autonomous systems. Another approach, the one that we develop 
here, is through geometric invariant manifold theory ([3, 8, 12]). There is a 
manifold M^, invariant for the full dynamics, which attracts all near-enough 
solutions, with an asymptotic phase property. The system restricted to the 
invariant manifold M^ is a regular perturbation of the fast (e=0) system. As 
remarked in Theorem 1.2 in Hirsch's early paper [4], a C^ regular perturba- 
tion of a flow with eventually positive derivatives also has generic convergence. 
So, solutions in the manifold will be generally well-behaved, and asymptotic 
phase implies that solutions track solutions in Mg, and hence also converge 
to equilibria if solutions on Mg do. A key technical detail is to establish that 
the tracking solutions also start from the "good" set of initial conditions, for 
generic solutions of the large system. 

For simplicity, we discuss here only the case of cooperative systems (niono- 
tonicity with respect to the main orthant), but proofs in the case of general 
cones are similar and will be discussed in a paper under preparation. 



2 Statement of Main Result 

We are interested in systems in singularly perturbed form: 

f = /(:^,y) (1) 

dy , . 
e^ = 5(-,y), 

where x £ M", y G R™, < e <C 1, and / and g are smooth functions. We will 
present some preliminary results in general, but for our main theorem we will 
restrict attention to the case when g has the special form g{x,y) = Ay + h{x), 
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where A is a Hurwitz matrix (all eigenvalues have negative real part) and h 
is a smooth function. That is, we will specialize to systems of the following 
form: 

§^ /(...) (2) 

^^ ^Ay + h{x). 
at 

(We remark later how our results may be extended to a broader class of 
systems.) Setting e to zero, we have: 

— ^ f{x,mn{x)), (3) 

where mo{x) = —A^^h{x). As usual in singular perturbation theory, our goal 
is to use properties of the limiting system (3) in order to derive conclusions 
about the full system (2) when < £ <IC 1. In this paper, A C B means that 
^ is a strict subset of B, while A C B contains the case oi A = B. We will 
assume given three sets K, K , and L which satisfy the following hypotheses 
(some technical terms are defined later): 

HI The set K is an n-dimensional C°° simply connected compact manifold 

with boundary. 
H2 The set L is a bounded open subset of K™ , and 

Mo = {{x, y) \y^ mo{x), x e K}, 

the graph of mo, is contained in K x L. 

H3 The fiow {tpt} of the limiting system (3) has eventually positive deriva- 
tives on K. 

H4 The set K is convex, and therefore it is p-convex too. 

H5 For each e > sufficiently small, the forward trajectory under (2) of 
each point in I? = IntK x L is precompact in D. 

H6 The equihbrium set Eo = {x G IntK \ f{x,mo{x)) = 0} is countable. 

H7 The set K C IntK is compact, and for each e > sufficiently small, the 
set D ~ K X L is positively invariant. 

Note that the equilibria of (2) do not depend on e, and the ones in D are 
in 1-1 correspondence with elements of Eo- The main theorem is: 



4 Liming Wang and Eduardo D. Sontag 

Theorem 1 Under assumptions H1-H7, there exists e* > such that for 
each < e < e* , the forward trajectory of (2) starting from almost every 
point in D converges to some equilibrium. 

Remark: A variant of this result is to assume that the reduced system (3) 
has a unique equilibrium. In this case, one may improve the conclusions of the 
theorem to global (not just generic) convergence, by appealing to results of 
Hirsch and others that apply when equilibria are unique. The proof is simpler 
in that case, since the foliation structure given by Fenichel's theory (see below) 
is not required. In the opposite direction, one could drop the assumption of 
countability and instead provide theorems on generic convergence to the set 
of equilibria, or even to equilibria if hyperbolicity conditions are satisfied, in 
the spirit of what is done in the theory of strongly monotone systems. D 



3 Terminology 

The following standard terminology is defined for a general ordinary differen- 
tial equation: 

where F : M.^ -^ M.^ is a C^ vector field. For any z G M^, we denote the 
maximally defined solution of (4) with initial condition z hy t —>■ (j}t{z),t £ 
/(z), where I{z) is an open interval in R that contains zero. For each t G M, 
the set of z G M^ for which 4)t{z) is defined is an open set W{t) C R^, and 
(j)t : W{t) -^ W{—t) is a diffeomorphism. The collection of maps {4)t]ti^R is 
called the flow of (4). We also write just z{t) for the solution of (4), if the initial 
condition z(0) is clear from the context. The forward trajectory of z G M^ is 
a parametrized curve t — » <f>t{z)- Its image is the forward orbit of z, denoted 
as 0+(z). The backward trajectory and the backward orbit 0-{z) are defined 
analogously. A set U C R is positively (respectively, negatively) invariant if 
0+{U) C U . It is invariant if it is both positively and negatively invariant. 

We borrow the notation from [8] for the forward evolution of a set U C 
F C R^ restricted to V: 

U -vt^ {Mp) ■ peU and (t>s{p) G F for all < s < t}. 

Let us denote the interior and the closure of a set U as IntC/ and U respectively. 
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Definition 1 The flow {04 } of (4) is said to have eventually positive deriva- 
tives on a set V C R^ if there exists ig such that -J^{z) > for all t > t^, 

zeV. 

If (4) is of dimension one, i.e. A^ = 1, {(j)t} has eventually positive derivatives 
automatically. In practice, the following sufficient condition is easier to check. 
If the vector field of (4) satisfies ^^{z) > 0, for all z e V, i =/= j, and the 

matrix ^{z) is irreducible for all z ^ V, then {4>t} has eventually positive 
derivatives. (This condition is not necessary.) 

Definition 2 An open set W C M^ is called p-convex if W contains the 
entire line segment joining x and y whenever x,y G V and x < y, where 
X < y means Xi < yi for all i ~ 1,- ■ ■ , N . 

The next lemma is a restatement of theorem 4.4 in [5]: 

Lemma 1. Suppose that the open set W C i?" is p-convex and the flow {4>t} 
of (4) has eventually positive derivatives on W . Let W^ Q W he the set of 
points whose forward orbit has compact closure inW. If the set of equilibrium 
points is countable, then z{t) converges to a equilibrium as t —>^ cx), for almost 
every z G W^. 

The following fact follows from differentiability of solutions with respect to 
"regular" perturbations in the dynamics; see [5], Theorem 1.2: 

Lemma 2. Assume V C W is a compact set in which the flow {</>*} has 
eventually positive derivatives. Then, there exists S > with the following 
property. Let {ipt\ denote the flow of a C^ vector field G such that the C^ 
norm of F{z) — G{z) is less than S for all z in V . Then there exists i* > 

such that if t > t^, and ips{z) G V for all s G [0,t], then -^{z) > 0. 

The Appendix reviews the definition of a "C (1 < t" < cxd) manifold M 
with boundary" in the sense used on geometric singular perturbation theory 
We denote the boundary of such a manifold M as 9M, and denote M\dM as 
IntM, when there is no confusion with the notation for the interior of a set. 

Definition 3 A compact, connected C" manifold M C M^ with boundary is 
said to be locally invariant under the flow of (1) if for each p in IntM there 
exists a time interval Ip = (^1,^2); for some ti < < t2, such that 4>t{p) G M 
for all t E Ip. 
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When e ^ 0, we can "stretch" the time (r = -), and consider the fast system: 

— ^sf{x,y) (5) 

The above system is equivalent to (1). The corresponding fast system for (2) 
is 

— ^ef{x,y) (6) 

-^ ^Ay + h{x). 

ClT 

Definition 4 Let M be an n- dimensional manifold (possibly with boundary) 
contained in {(x,y) | g(x,y) — 0}. We say that M is normally hyperbolic 
relative to (5) if all eigenvalues of the matrix -^ (p) have nonzero real part for 
every p £ M . 



4 Proof of the Main Theorem 

Recall the definition of Mq — {{x,y) \ y = mo{x),x G K}. Since K is an 
n-dimcnsional C°° compact manifold with boundary, and ttiq is a smooth 
function, Mq is also an n-dimensional C'°° compact manifold with boundary. 
Our proofs are based on Fenichel's theorems [3], in the forms presented 
and developed by Jones in [8]. 

Fenichel's First Theorem Under assumption HI, if Mq is normally hyper- 
bolic relative to (6), then there exists eg > 0, such that for every < £ < £o 
and r > 0, there is a function y — m^{x), defined on K , of class C"" jointly in 
X and e, such that 

Me = {{x, y) I y = TOe(a;), x <^ K] 

is locally invariant under (2), see Figure 1. 

The requirement that Mq be normally hyperbolic is satisfied in our case, 
as g{x,y) = Ay + h{x) and therefore -^{p) = A, which is Hurwitz, for each 
p€Mo. 
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K ^ 

Fig. 1. For simplicity, we sketch manifolds AL and Mo of a system where n — m = 1. 
The set K is & compact set in x, and L is an open set in y. The red curve denotes 
the locally invariant manifold Me and the black curve denotes Mo. 



We will pick a particular r > 1 in the above theorem from now on. 

Let us interpret local invariance in terms of equations. Let {x(t),y{t)) be 
the solution to (2) with initial condition {xo,yo), such that xq £ IntK and 
2/0 — 'rng{xo). Local invariance implies that {x{t),y{t)) satisfies 



dx{t) 
dt 



f{x{t),m,{x{t))) 



y{t) = me{x{t)), 



(7) 

(8) 



for all t small enough. Actually, this is also true for alH > 0. The argument is 
as follows. By H5, (x{t), y{t)) is well-defined and remains in D for all t > 0. 
Let T = {i > I y{t) = ms{x{t))}. Then, T is not empty, and T is closed 
by the continuity of m^{x{t)) and y{t). Also, T is open, since M^ is locally 
invariant. So T = {i > 0}, that is, x{t) is a solution to (7) and y{t) — m^{x{t)) 
for all t > 0. 

In (7), the x-equation is decoupled from the y-equation, which allows us 
to reduce to studying a lower-dimension system. Another advantage is that, 
as e approaches zero, the limit of system (7) is system (3), which describes 
the flows on Mq. If Mq has some desirable property, it is natural to expect 
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this property is inherited by the perturbed manifold Mg. An example of this 
principle is provided by the following lemma. 

Lemma 3. Under assumptions H1-H3, for each < e < Sa, the flow {"0*} 
of (7) has eventually positive derivatives on IntK . 

Proof. Applying Lemma 2, there exist 5 > Q such that when the C^ norm 
of mo(x) — rn^{x) is less than 5 for all x & K, there exists i, > with the 

property that: if t > t, and ipsix) G K for all s G [0,t], then -^{x) > 0. Since 
rUe is of class C", jointly in x and e, we can pick e > small enough to control 
5. If we then prove IntJC is invariant under (7), that is, for any xo G Intii', the 
solution x{t) of (7) with initial condition xq stays in hiiK for all t > Q, then 
we are done. Let us now prove that IntK is positively invariant under (7). 

Let yo = fTisixo) a-iid y{t) = me{x{t)). Then, {x{t),y{t)) is the solution to 
(2) with initial condition (a;o,J/o) G D. By H5, {x(t),y{t)) stays in D for all 
i > 0, and therefore x{t) G Intii' for aU t > 0. D 

Flows with eventually positive derivatives have particularly appealing 
properties, as in Lemma 1. To apply that lemma, we need to check two con- 
ditions. First, for every point in IntK, its forward trajectory under (7) has 
compact closure in IntK. Second, the number of equilibria of (7) is count- 
able. Suppose that the first property does not hold, and let x{t) be a solution 
to (7) with x(0) G IntK but \inij ^ oc x{tj) ^ IntK for some sequence {tj}. 
So, yx{t),me{x{t))) is a solution for (2), and its forward orbit is not precom- 
pact in D. This violates H5. To check the second condition, wc introduce the 
following sets: 

E,^{xe IntK I f{x,m,{x)) = 0} , < £ < 1. 

We claim that Eg C Eq for all e small enough, which implies that E^ is 
countable, by H6. Let us prove the claim. It is clear that Eq consists of the x- 
coordinates of all equilibria of (2) in D. Fix < e < £o and pick any xq <E E^, 
2/0 = iTieixo)- The solution (x(t),j/(t)) to (2) with initial condition (xo,yo) 
satisfies (7) and (8) for t small enough. But 

^ = f{x{t),y{t)) = /(x(i),m,(x(t))) EE 0, 

so (x{t),y(t)) ^ {xo,yo) for all t > 0, and therefore xq G Eq- Applying Lemma 
1 we have: 
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Lemma 4. Under assumptions H1-H6, for each < e < Sq, there exists a 
set C^ C IntK such that the forward trajectory of (7) for every point of C^ 
converges to some equilibrium, and the measure of IntK \ C^ is zero. 

Until now, we have discussed the flow only when restricted to the locally 
invariant manifold M^. The next theorem, stated in the form given by [8], 
deals with more global behavior. In [8], the theorem is stated for e > 0, but 
some properties also hold for e = ([9]). (We will apply this result again with 
a fixed r > 1.) The notation [—5, S\ stands for the cube { (j/i, . . . , y™) | — 5 < 

\yi\<5}. 

Fenichel's Third Theorem Let £o be as in Fenichel's First Theorem. Under 
assumption HI, if Mq is normally hyperbolic relative to (6), then there exists 
< ei < £o 0,1^^ 5 > Q such that for every < £ < ei and r > 0, there is a 
function 

h^: K X [-S, S] -^ M" 

such that the following properties hold: 

1. For each x G K, he{x,0) — x. 

2. The image of the map 

T^:Kx [-6, (5] ^ M" X R" 

(x. A) H^ {hs{x, A), A + m^{h^{x, A))) 

is defined as the stable manifold W^{M^) of M^. For p = (x,TOe(x)) G M,,, 
the stable fibers W^{p), defined as T^dx} x [—(5,(5]), form a "positively 
invariant" family when e ^ Q, in the sense that 

W:{p)-WS^Me)t^W:{ct^t{p)). 

3. "Asymptotic Phase". There are positive constants k and a such that for 
any p,qe M"+™, if q G W^{p), e^O, then 

\Mp) ~ Mq)\ < ke-^' 

for all t > as long as 4>t(p) and 4>t((l) stay in W^{Me). 
4- The stable fibers are disjoint, i.e., for qi G W^{pi), i = 1,2, either 

W^ipi) n W/(p2) = or W/(pi) = Wiip2). 
5. The function h^{x, X) is C^ jointly in e, X and X. Whene ~ 0, ho^six, X) = x. 
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W^((]) ,(p)) 




Fig. 2. To illustrate the geometric meaning ol Fenichel's Third Theorem, we sketch 
the locally invariant manifold and stable fibers of a system, in the case of n=m=l. 
The dimensions of the manifolds M^, Mo, and stable fibers are one. M^ is the graph 
of A = 0, and Mo is the graph of mo (x) — rrie {x) (black curve) . These manifolds may 
intersect at some equilibrium points. Through each point p G M^ (x-axis), there is 
a stable fiber W" (p) (blue curve) . We call p the "base point" of the fiber. The fiber 
consists of the pairs {he{x, A), A), where [A] < (5. If a solution (purple dashed curve) 
starts on fiber W^(p), after a small time t, it evolves to a point on another stable 
fiber W"{cf>t{p)) (light blue curve); this is the "positive invariance" property. 



The next lemma gives a sufficient condition to guarantee that a point is 
on some fiber. 

Lemma 5. Let ei and S be as in Fenichel's Third Theorem. There exists 
< £2 < £i, such that for every < s < 62, the set 



As := {{x, y) \ X e K, \y - mo{x)\ < -} 

is a subset ofW^{Me)- 

To prove this lemma, we need the fohowing result: 

Lemma 6. Let U and V be compact, convex subsets of M" and 
tively. Suppose given a continous function 



respec- 
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(/. : U xV ^W xW"- 

{x,y) t-> {(pi{x,y),(l)2{x,y)) 

satisfying \\(j)i{x,y) - x\\ < pi, \\(j)2{x,y) - y\\ < p2 for some pi > 0, p2 > 
and all {x,y) G U xV . Then every point (a, (3) (z U xV with dist{a, dU) > pi 
and dist{j3, dV) > p2 is in the image of (f>. 

Proof. For such a point (a, /3), consider the map <?(a;, y) = {(Pi(x,y),<p2{x,y)) := 
{x,y) — {(l)i{x,y),(l)2{x,y)) + {a,j3). Thus 'P maps U xV into itself. If not, 
say <I>i{x,y) is not in U, that is, x — (j)i{x,y) + a is not in U. Since 
||a; — (j)i{x,y)\\ < pi, then d[st{a,dU) < pi, contradiction. The case when 
$2{x,y) is not in V follows similarly. Since <? maps U x V into itself, and 
the product of convex sets is still convex, by Brouwer's Fixed Point Theo- 
rem, there is some {x,y) dzU xV so that <P{x,y) ~ {x,y), which means that 
{(f>i{x,y), (f>2{x,y)) = (a,/3), as we wanted to prove. D 
Proof of Lemma 5. Define for each < e < ei, the map 

(j)e- Kx [-(5, 5] -^ K" X K" 

{x, A) 1-^ {he{x, A), A + m^{he{x, A)) - mo(/ie(a;. A))). 

By property 5 in Fenichel's Third Theorem and the compactness of K, we 
have 

\\h,{{x,m,{x)),\)-x\\<Ci(e), 

and 

|1[A + m,{h,{x, A)) - m„{K{x, A))] - A|| < C2{e) 

for some positive functions Ci of e such that Ci ^ as e ^ 0, i = 1, 2. For any 
such £, apply Lemma 6 with U ^ K,V ~ [—5, 5], px ~ Ci(£), p2 ~ C'2(£) and 
(f) = 4>e- Since dist(iir, dK) and | are independent of e, we can pick < £2 < £i 
such that dist(if , dK) > Ci(e) and | > C2(e) for all e G (0, £2)- By Lemma 6, 

Kx[-^-,^-]cMKx[-S,S]) (9) 



Define a map 



TTo : IR" X K'" -^ K" X M'" 

(a;. A) ^ {x,X + mo{x)) 
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Consider the composition of ttq and 0^: By property 2 in Fenichcl's Third 
Theorem, its image, ttq o (I)^{K x [—(5,(5]), is the stable manifold W^{Mg) of 
Mg. According to (9), ttq{K x [-^,^]) C W^iM^). Notice that As ^ TToiK x 
[— |, |]), we are done. D 

Lemma 7. Let e-i he as in Lemma 5, and 5 as in Fenichel's Third Theorem. 
Under assumption H7, there exists < £3 < £2 such that for each < £ < £3, 
if p & D, then there exists Tq > 0, and 4>t{p) G As for all t >T(). 

Proof. Setting z = y — too(x) and r = -, (2) becomes 

-— = sf{x,z + mQ{x)) 
dr 

dz 

— = Az — emL(x)f(x, z + mo(x)). 

dr 

So 

z(t) = 2(0)6-^^ -e f e'^^''-'''>mQ{x)f{x, z + mo(x)) ds. 
Jo 

Notice that l|e"^*|| < Ce'^*, for some positive constant C and negative /3, which 
is greater than the real part of all eigenvalues of A. So, 



£ / e^^'^ ^^m'f^{x)f(x,z + mo{x)) ds 



2eMC 
/o 

where M is an upper bound of the function ||?7iQ(x)/(x,y)|| on D. Let 



"^' and T-max{4lnl^M,0}. 



SMC 



Then, we have ||z(t)|| < | for all r > Tq. Back to the slow time scale, we let 
To == eT(^. Therefore, (j)t{p) £ As for all t > To, derived from H7. 

Remark: Except for the normal hypcrbolicity assumption. Lemma 7 is the 
only place where the special structure (2) was used. Consider a more general 
system as in (1), and assume that g{x,mQ{x)) = on K for some smooth 
function toq. By the same change of variables as in the above proof, (1) is 
equivalent to 
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— = ef{x,z + mo{x)) 

dz 

-— — g(x,z + mQ(x)) - em'Jx)f{x, z + mo(a;)). 

CLT 

The only property that we need in the lemma is that for any initial condition 
(a:;(0), z(0)), the solution {x(t),z{t)) satisfies 



limsup|z(t)| < 7 limsuprf(f) 

t — *oo \ t — *oo 

where 7 is a function of class /C, that is to say, a continuous function [0, cx)) —^ 
[O,cxo) with 7(0) = 0, and d{t) = eTOg(a;(i))/(x(i), z(t) +mo(a;(t))). In terms 
of the functions toq and g, we may introduce the control system dz/dt — 
G{d(t), z) + u{t), where d is a compact-valued "disturbance" function and u is 
an input, and G{d, z) = g(d, z + TOo((i)). Then, the property of input-to-state 
stability with input u (uniformly on d) , which can be characterized in several 
different manners, including by means of Lyapunov functions, provides the 
desired condition. D 

Lemma 7 proves that every trajectory in D is attracted to As and therefore 
is also attracted to M^. This will lead to our proof of the main theorem. 

Proof of the main theorem. 

Choose £* — £3, defined in Lemma 7. For any p ^ D, there arc three cases: 

1. p £ Mg. By Lemma 4, the forward trajectory converges to an equilibrium 
except for a set of measure zero. 

2. p £ As C W^{Ms). Then p is on some fiber, say W^{p), where p — 
{x,m^{x)) G Mg. If X is in Cg (defined in Lemma 4), then (j)t{p) —>■ q, 
for some q E Eq. By the "asymptotic phase" property of Fenichel's Third 
Theorem, (jjtip) also converges to q. To deal with the case when x ^ C^, it 
is enough to show that the set 

xelntK\Ce 
as a subset of M'"+" has measure zero. Define 

Te= flnt^\Ce) X [-6,5]. 
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Since IntK\C^ has measure zero in R", also J> has measure zero. On the 
other hand T^{J-'g) — B^, and Lipschitz maps send measure zero sets to 
measure zero sets, we are done. 
3. p e D \ ^5. By Lemma 7, (/)t(p) G As for all t > Tq. Without loss of 
generality, we assume that Tg is an integer. If (f)T„ (p) G AsXBg, then (f)t (p) 
converges to an equilibrium. Otherwise, p £ lJfc>o fcez 4'-k(Be)- Since the 
set Be has measure zero and (jj-k is Lipschitz, (j)-k{Be) has measure zero 
for all fc, and the countable union of them still has measure zero. D 



5 An Example 

Consider the following system: 

-^ ^lt{yi,---,yrn) - (3i{xi,...,Xn), i=l,...,n, 

e-j^ = -djyj -aj{xi,...,Xn), dj > 0, j==l,...,m, (10) 

where aj, [3i and 7^ are smooth functions. We assume that 

1. The reduced system 

-^ =Ji{- — ,...,--^)-f3^{xi,...,Xri) -.^ F^{xi,...,Xn), i^l,...,n 
has partial derivatives that satisfy: 

— = V- — — — - — > Ofor i^^ fc (11) 

dxk ^ di dyi dxk dxk 

2. For each i, 

lim min (3i{x) — +00 (12) 

«^+oc rreSj(u) 



and 



lim max (3i{x) = — cxo (13) 



where Si{u) is the set of vectors in R" whose ith coordinate is u. (For 
n = 1, this means simply that lima;^±oo (3i{x) = ±00.) 
3. There exists a positive contant Mj such that \aj{x)\ < Mj for all x £ R". 
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4. The number of roots of the system of equations 

7i(ai(x), . . . ,a„i{x)) = (ii{x), i = 1, . . . ,m 

is countable. 

We are going to show that on any large enough region, and provided that 
£ is sufficiently small, almost every trajectory converges to an equilibrium. To 
emphasize the need for small e, we also show that when e > 1, a limit cycle 
could appear. 

To apply our main theorem, we take 

L = {yeM" I \y,\<b,, j = l,...,m}, 

where bj is an arbitrary positive number greater than -r^. Picking such bj 

assures 2/j-J^ < for all a; G M and \yj\ = bj, i.e. the vector field points 
transverally inside on the boundary of L. Let 

K — {x eW^ \ - ai^2 < Xi < ais, i = 1, . . . , n } 

where a^.i and 0^.2 can be any positive numbers such that 

f3i{x)>Ni:^ max |7j(j/i, . . . ,y„) 

\yi\<''i 

whenever x G M" satisfies that its ith coordinate Xi > Ui^i, and 

whenever x G K" satisfies that its ith coordinate Xi < —ai^2- All large enough 
Oi.j's satisfy this condition, because of the assumption made on /3. 

So, we have Xi-j^ < for all y G L, Xi — a^^i and Xi — —ai^2- We then 
take ^ 

K = {x eM."^ \ - ai^2 - i < Xi < fli^i + 1, 2 = 1, . . . n }, 

D ~ KxL and D — IntK x L. Thus, the vector field will point into the interior 
of D and D. Hypotheses H5 and H7 follow directly from this fact. (Sketch: 
H7 is obvious. Suppose H5 does not hold. Then, there exists some solution 
{x{t),y{t)) of (10) in D, and a sequence {tj} such that {x{tj),y{tj)) -^ {x,y) 
as j -^ cxD. Suppose that yk — bk for some k G { 1, . . . ,?n}, and {yk{tj)} is 
strictly increasing to bu- This will contradict the fact that -^ < above the 



16 Liming Wang and Eduardo D. Sontag 

t/fc-nuUclinc. The other cases foUow similarly.) H3 follows from our assumption 
1, and it is easy to see the other hypotheses also hold. By our main theorem, 
for sufficiently small e, the forward trajectory of (10) starting from almost 
every point in D converges to some equilibrium. 

On the other hand, convergence does not hold for large e. Let 

n = l, (3i{xi) = -^ - xi, m == 1, ai{xi) = 4tanha;i, ^{yi) = yi, di ^ 1. 

It is easy to verify that (0,0) is the only equilibrium. When £ > 1, the trace 
of the Jacobian at (0, 0) is 1 — i > 0, its determinant is — > 0, so the 
(only) equilibrium in D is repelling. By the Poincare-Bendixson Theorem, 
there exists a limit cycle in D. 

Remark: The conditions (11), (12), and (13) are satisfied, in particular, if 
one assumes the following easier to check conditions on the functions /3i's, 
a^'s, and 7i's, The functions (3i are asked to be so that: 

oxk 

for every i,k — 1, . . . , n such that i ^ k (cooperativity condition among Xi 
variables), and also so that: 

lim Pi{xi,. . . ,Xn) = +00 (14) 

Xi — ^-\-00,....Xn — * + 00 

lim (3i{xi,. . . ,Xn) ^ -00. (15) 

xi — > — 00, ...,3:^1 — * — 00 

These last conditions are very natural. They are satisfied, for example, if 
there is a linear decay term —XiXi in the differential equation for each Xi, 
and all other variables appear saturated in this rate. Since -q^{x) < for all 
i T^ k, (14)-(15) imply that conditions (12) and (13) both hold. Regarding the 
remaining functions, we ask: 



E 



djt dai _^ p 



,^^ dyi dxk 

for alH, fc = 1, . . . , n such that i ^ k. This condition can be guaranteed to hold 
based only upon the signs of the partial derivatives: it holds true if there is no 
indirect negative effect (through the variables yi) of any variable Xk on any 
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yi - X2 

Fig. 3. Example 

other variable xi. The diagram shown in Figure 3 illustrates one such influence 
graph (signs indicate signs of partial derivatives), for n = m = 2. Observe that 
this example cannot describe a monotone system (with respect to any orthant 
cone, i.e., it is not cooperative under any possible change of coordinates of 
the type Xi —> —Xi or yi —^ —yi). An entirely analogous example can be done 
for any n — m, the key property being that each variable Xi "represses" its 
associated variables yi and the j/j's "enhance" some or all other variables. 



6 Appendix 

Definition 5 The closed half-space iJ' C M', is defined as follows: 

H^ = {(xi,a;2,--- ,xi) £ R' | xi > 0}. 
The boundary of H\ denoted as dH^ , is R'^^. 

Definition 6 A subset M C M is called a l-dimensional C" manifold with 
boundary if it satisfies: 

1. There exists a countable collection of open sets V"' C M. , a ^ T, where 
T is some countable index set, so that, with U" = V"' P| M , one has M — 

2. There exists a C^ diffeomorphism x" defined on each U°' which maps C/" 
to some set Wf^H^ where W is some open set in M'. 

The boundary of M , denoted as dM , is the set of points in M that are mapped 
to dH^ under x", for some a G Z. 
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A compact manifold is a manifold that is compact as a topological space. The 
definition implies that dM is a well defined C manifold of dimension I — \, 
and M \ dM is an ^dimensional C manifold; see [10] for the details. 
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